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Abstract. We study coarsening phenomena in three different simple exclusion 
processes with quenched disordered jump rates. In the case of the totally asymmetric 
process, an earlier phenomenological description is improved, yielding for the time 
dependence of the length scale ~ t/(\nt) 2 , which is found to be in agreement 
with results of Monte Carlo simulations. For the partially asymmetric process, the 
logarithmically slow coarsening predicted by a phenomenological theory is confirmed 
by Monte Carlo simulations and numerical mean-field calculations. Finally, coarsening 
in a bidirectional, two-lane model with random lane-change rates is studied. Here, 
Monte Carlo simulations indicate an unusual dependence of the dynamical exponent 
on the density of particles. 
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1. Introduction 

Single-file transport is known to be susceptible to heterogeneities of the track on which 
the process takes place — one may think of formation of jams in vehicular traffic p~| 
or shocks in the intracellular traffic of molecular motors on the filamentary network 
[2]. Such systems are frequently modeled by some variant of the asymmetric simple 
exclusion process (ASEP) [31 0], where particles jump stochastically on a lattice and 
interact through excluding multiple occupations on lattice sites. For this model, in 
its homogeneous form, many exact results are available [51 Ej and, therefore, it has 
become a paradigmatic model of non-equilibrium systems. Motivated by understanding 
transport in heterogeneous systems, the ASEP with quenched random jump rates has 
been the subject of ongoing research [Zl El El HDl IHl H21 HS1 H31 HS1 HSl IIZI HB] - Besides 
the experimental relevance, the effect of disorder on this model is also a challenging 
problem as general exact solutions are still lacking here, and the bulk of the results 
have been obtained by mean-field approximations, phenomenological theories based on 
extreme- value statistics and Monte Carlo simulations. 

Since the ASEP can be mapped to a growth process of a 1 + 1 dimensional 
surface [19] which is described by the Kardar-Parisi- Zhang equation [20J, results for 
exclusion processes are also important to understand the fundamentals of surface physics 
[2Tj or directed polymers in random media [20]. Due to the steady current, the 
correlation length diverges in these systems and, consequently, they exhibit a non- 
equilibrium critical scaling behavior that can be classified into universality classes [22]. 
Understanding the effects of disorder on lattice gas models thus reveals the behavior 
of the corresponding interfaces. In particular, variants of the ASEP with quenched 
randomness in the hopping rates are related to surface growth models with columnar 
disorder (see [22]). 

Typically, these models relax toward the steady state much slower than the 
corresponding homogeneous ones; in certain cases, the coarsening dynamics are 
logarithmically slow. Owing to this, confirming the predictions of phenomenological 
theories by Monte Carlo simulations is a very hard numerical task. As it has been 
recently shown, inhomogeneous exclusion processes can be parallelized and simulated 
efficiently on a GPU [23], similar to other interacting systems [211 1251 [26] . In this 
work, we shall get use of this new technique (among other methods, such mean- 
field approximation) in studying the dynamics of three variants of the disordered 
ASEP: (i) the disordered totally asymmetric simple exclusion process (DTASEP), where 
particles move unidirectionally, (ii) the disordered partially asymmetric simple exclusion 
process (DPASEP), where particles can move in both directions and (iii) a recently 
introduced bidirectional two-lane model, where the disorder appears in the lane-changes 
[27] . We shall see that the above disordered processes with parallel (or synchronous) 
update procedure have essentially the same physical properties as those with (the 
most frequently studied) random sequential (or asynchronous) update procedure. In 
comparing the results with the mean-field description we shall enlighten the limitations 
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of the latter in case slow modes related to moving domain walls are present. 

The paper is organized as follows. In section [2] and [3], the coarsening phenomenon 
in the DTASEP and in the DPASEP is investigated, respectively. Section H] is devoted 
to the steady state and the coarsening in the bidirectional two-lane model. Finally, 
results are discussed in section |5j Some technical details concerning the parallelization 
are given in the Appendix. 

2. Coarsening in the disordered totally asymmetric simple exclusion 
process 

First, we shall study the disordered totally asymmetric simple exclusion process defined 
as follows. A one- dimensional periodic lattice with L sites is given, on which N particles 
reside, at most one on each site. The TASEP is a continuous-time stochastic process 
in which particles move from site % to site % + 1 independently with rate pi provided 
site i + 1 is empty. The jump rates Pi are independent random variables drawn from a 
bimodal distribution: 

f{pi) = {l-<j>)S[p i -l) + <f>S{p i -r) , (1) 

where < r < 1 and < < 1. Although exact solutions are lacking for this model, 
the main features of the steady state and the non-stationary (coarsening) behavior are 
well understood by means of a phenomenological description [121 HSj- This rests on 
that particles accumulate behind consecutive stretches of bonds with rate Pi = r (called 
bottlenecks) and form a high density domain where the local density is greater than 1/2. 
Provided that the global density N/L is close to 1/2, and the system is initiated in a state 
with a uniform density, then high-density segments form behind bottlenecks and start to 
grow as time elapses. During this process, segments at short bottlenecks lose particles in 
favor of those at longer ones until the former gradually vanish. Ultimately, provided the 
system was finite, a single, macroscopic high density domain accumulates behind the 
longest bottleneck while in the rest of the system the local density is below 1/2. So, the 
system undergoes a coarsening process which is characterized by the dependence of the 
typical size of segments £(t) on time. In a simplified description of this phenomenon, the 
bottlenecks are treated to be independent and are characterized by a (length dependent) 
maximal current. Furthermore, the essence of the coarsening process is captured if 
one concentrates merely on the evolution of length of high density segments behind 
bottlenecks. Under these assumptions, the model reduces to a disordered zero-range 
process [28], where the dynamical exponent z describing the coarsening as 

at) ~ (2) 

is known to depend on the asymptotics p(w) ~ (w — c) n of the distribution of jump 
rates at the lower edge c as z = (n + 2)/(n + 1) [29]. In the DTASEP, the maximal 
current through a bottleneck of length I is [30] 

J(l) =r/4 + 0(r 1 ) (3) 
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and the distribution of the bottleneck length is geometric. This yields formally n = oo 
and z — 1. In this special limit, non-rigorous arguments based on extreme value statistics 
indicate that the coarsening is not simply linear but logarithmic corrections arise as 
f(t)~t/lnt ETJ. 

Here, we shall reconsider this reasoning and refine it by applying the statistics of 
extremes in a more precise way. Let us consider a growing high-density segment of length 
Clearly, the rate of growth of £j(t) is determined the difference AJj between the 
current Jj_i through the bottleneck on the left hand side of the growing segment and 
the current Jj through the bottleneck on its right hand side: 

ii(t) ~ AJ t . (4) 

The two bottlenecks are separated by a distance 0(£(t)) and, in this domain, the 
bottleneck on the right hand side is the longest one (having length and the other one 
is the second longest one (having length I2). Using Eq. ([3]), the current difference can 
be written in terms of bottleneck lengths as 

&J i ~l/l 2 -Vh = {h-h)lhh. (5) 

Approximating, for the sake of simplicity, the geometric distribution of bottleneck 
lengths by the continuous exponential distribution P<(1) = 1 — e~ la ^^^ 1 , it is 
straightforward to show that the most probable values of l\ and 1% among £ events, 
l\ ~ ln£/ln(l/<2) and l\ ~ ln(£/2)/ln(l/Q), respectively, are shifted by a finite value. 
Furthermore, the variances of the distributions of l\ and l 2 tend to finite values in the 
limit £ 00. Therefore the limit distribution of l± — I2 also has a finite most probable 
value and variance and, consequently, the typical current difference scales with £ as 

AJ<~ l/Z^^l/Cln 2 (6) 
for large £. Putting this into Eq. (T4J) and integrating yields for the growth of £(t) 

m-t/(\n t y (?) 

in leading order for long times. 

For measuring the typical length of segments we have made use of a well-known 
relationship between the ASEP and a simple one-dimensional surface growth model 
where height differences A/i, = h i+ i — hi are related to the occupation number rii of the 
ASEP as Ahi = 2n^ — 2N/L [19]. A usual measure of the width of the surface is given 
by 

w{^t) = [\j:h>{t)-(\j:h l {t)f] 1,2 ■ ( § ) 
l 1=1 l i=i 

In case of the coarsening DTASEP, the equivalent surface consists of roughly linearly 
ascending and descending parts corresponding to high and low density segments, 
respectively. Therefore the quantity given in Eq. (JS} is proportional to the typical 
value of 

We performed numerical simulations of the model both with random sequential 
updates and parallel sub-lattice updates (see the Appendix) and used the bimodal 
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Figure 1. Dependence of the surface width of the DTASEP on time, obtained by 
Monte Carlo simulations with parallel updates for L = 16000, 48000. The number of 
independent samples are given in parentheses. 




Figure 2. Scaling plot of the width W(L,t) obtained for different system sizes by 
Monte Carlo simulations with parallel updates. 



distribution with parameters <fi = 1/2 and r = 1/4. The density of particles was 
N/L = 1/2. We have considered systems with size L = 16000 and L = 48000 and 
measured the time evolution of the surface width. The average of W(L,t) is plotted 
against t/(\nt) 2 in Fig. CD Here, a linear dependence can be seen for not too long 
times, where finite size effects are negligible. Data for finite sizes satisfactorily follow 
the scaling law 

W(L,t) = LW{t/{\ntf/L), (9) 

see Fig. [2j although, with strong corrections to scaling. 

In order to compare the prediction in Eq. (ED) with numerical results more precisely 
we have calculated the effective dynamical exponent by l/z e g(t) = dba.W(t,L)/d\nt. 
Due to the logarithmic factor in Eq. (j7|), the leading order correction term in the 
finite-time effective exponent is slowly decaying, having the form 

l/2 eff (t) = l-2/lnt + --- (10) 
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Figure 3. Effective dynamical exponents plotted against time obtained by Monte 
Carlo simulations with random sequential updates (+,x) and by numerically solving 
the numerical mean field equations (*). The curves correspond to the predictions of 
the phcnomenological theory. 



in an infinite system. As can be seen in Fig. [3j for moderate time scales the effective 
exponents still deviate from the predictions of the simple phenomenological theory, 
nevertheless, the agreement with Eq. (fTUl) is satisfactory for long times. Due to the 
slow, logarithmic convergence, the estimates on the exponents obtained on small scales 
differ significantly from the asymptotic values, see the numerical results in Ref. [3Tj . 

Note, that the simplest scaling hypothesis for roughening surfaces, which is valid 
for, among others, the above model in the absence of disorder is of the form 

W{L,t) = L a W{t/L z ) (11) 

rather than that given Eq. ([9]). Here, according to the usual terminology, a is the 
roughness exponent characterizing the finite-size scaling of the stationary width of 
the surface, i.e. W(L,t — > oo) ~ L a , while the surface growth exponent (3 = zja 
describes the evolution of the surface width in an infinite system for long times as 
W(L — > oo,t) ~ t° . Hence, apart from the logarithmic correction, the surface growth 
model corresponding to the DTASEP is characterized by the exponents a — /3 — z — 1. 

We have also studied the model within a mean-field approximation in which the 
expected values (ni(t)n i+ i(t)) in the master equation of the process are replaced by 
{rii(t)) {n i+ i(t)) . Then the master equation turns into the following evolution equations 
for the local densities pi(t) = (ni(t)) 

^ = p i _ip i _i(l - pi) -piPi(l - pi+i), % = 1,2, . . . ,L. (12) 
dt 

This approximation gives the general features of the ASEP correctly, but, since the 
correlations are neglected, it may yield incorrect critical exponents. From the point of 
view of the above phenomenological theory an important difference is that the finite-size 
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corrections of the current in the mean-field model are in the order of l~ 2 |31j . rather 
than see Eq. ([3]). Repeating the above calculations otherwise unchanged, leads 
to a different outcome for the form of the logarithmic correction; one obtains namely 
~ t/(lnt) 4 and, correspondingly, l/z cS (t) = 1 — 4/ lnt + • • •. We have solved Eqs. 
(112j) numerically by using the 4th order Runge-Kutta method for 10 2 samples with 
L = 10 6 . The calculated effective exponents are compared with the predictions of the 
theory in Fig. [3j As can be seen, on the time-scales of the numerical investigation a 
deviation from the predictions is present, which is, however, decreasing with increasing 
time. 

3. Coarsening in the disordered partially asymmetric exclusion process 

The dynamics are much different in the quenched disordered partially asymmetric 
exclusion process (DPASEP), where particles can move in both directions with random 
rates Pi and qi to the right and to the left, respectively, from site i 0, [8j [HI EES]. Here, 
due to the large fluctuations of the random potential landscape [32], which is itself a 
random walk, the displacement of a single particle increases ultra-slowly as 

WW)~{\ntf^ (13) 

with the barrier exponent ip = 1/2, when the average force F = \n(pi/q i+ i) acting 
on the particle is zero [33]. Here, (•) denotes an average over different stochastic 
histories, whereas the overbar denotes an average over the random transition rates. 
In the following, we shall restrict ourselves to the unbiased case, i.e. F — 0. In the 
presence of many particles (i.e. in the DPASEP), numerical simulations showed that the 
motion of a tagged particle follows the same law given in Eq. ( fl3l) [8] and, in accordance 
with the activated dynamics, the stationary current in finite rings of size L behaves as 

- In | J(L) | ~ L*. (14) 

When the system starts from a homogeneous state it shows a coarsening phenomenon 
similar to the DTASEP. Queues form at potential barriers, where the local density is 
close to one, leaving the rest of the lattice almost completely empty. As time elapses, 
queues at smaller barriers gradually dissolve and particles accumulate behind larger and 
larger barriers. Scaling considerations based on the finite-size scaling behavior of the 
current in Eq. ( !T4|) lead to that the typical size of high-density segments increases with 
time in an infinite system as [H] 

£(t)~ln^[t/ln(t)]. (15) 

This means that the corresponding surface is characterized formally by the exponents 
a — 1, /? = and z = oo. We mention that, in the presence of a bias (F > 0), the 
scaling is normal, i.e. has the form given in Eq. [TTJ with a = 1 and with /3 and z 
continuously varying with F [T4"] . 




Figure 4. Dependence of the surface width of the DPASEP on time, obtained by 
Monte Carlo simulations for sizes L = 512,2048,4096,8192. 



Our aim here is to test the coarsening law in Eq. ( Ti~5]) by measuring the dependence 
of the corresponding surface width W{t) in Monte Carlo simulations. Numerical results 
obtained with the parallel updating procedure for different system sizes are shown in 
Fig. HI As can be seen, not considering long times where the data are affected by the 
finite size of the system, the results are compatible with the law given in Eq. ( fl5l) . We 
have also considered the model within the mean-field approximation, which has recently 
been demonstrated to describe both the stationary behavior and the dynamics of the 
model correctly [18J. The dynamical mean-field equations read for this model as 

= {Pi-iPi-i + Qipi+i)^ - Pi) ~ [Pi(l - Pi+i) + - Pi-i)]pi- (16) 

In the earlier work, the coarsening length scale has been studied within the mean-field 
approximation indirectly by measuring the typical distance between adjacent peaks in 
the density profile [18]. Here, we have solved the equations Eqs. fTTol) numerically and 
calculated directly the surface width W(t) as a function of time. Results shown in Fig. 
[5] are in good agreement with Eq. (fT5|) . 

4. Dynamics of a bidirectional, two-lane exclusion process with random 
lane-change rates 

4-1. Definition of the model and preliminaries 

In the final part of this work, we shall study a recently introduced bidirectional, two- 
lane exclusion process, where the lane-change rates are random variables [3U EZj- The 
precise definition of the model is given as follows. Two parallel, one- dimensional, 
periodic lattices are given, each of them with L sites. Particles move within a lane 
unidirectionally, following the dynamical rules of the TASEP with unit jump rates but 
the directions of motion in the two lanes are opposite. In addition two this, particles 
change lanes from site i of lane A(B) to the neighboring site of lane B(A) with rate Ui(vi) 
provided the target site is empty. The model and its possible transitions are illustrated 
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Figure 5. The same plot as in Fig. [4j but the data are obtained by solving the 
dynamical mean-field equations. The number of samples was 10 2 . 



in Fig. [6j Notice that the direction of motion of particles can be regulated in this model 



r\ r r 
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Figure 6. Sketch of the bidirectional, two-lane model under study. Allowed 
transitions with the corresponding rates are indicated. 



by tuning the asymmetry in the lane-change rates. This arrangement was intended to 
model traffic of molecular motors moving along oppositely oriented filaments, which can 
be realized in experiments and may be relevant with respect to in vivo processes [35] . 
A single motor in this environment with homogeneous and symmetric lane-change rates 
performs normal diffusion, however, with an enhanced diffusion coefficient compared to 
that of the symmetric random walk [36]. Recently, the model has been reinterpreted as 
a model of motors able to move bidirectionally along the filament |37j . 

We shall consider the disordered version of the above model where the pairs of rates 
(ui,Vi) are i.i.d. random variables. In case there is only a single particle in the system, 
one can introduce an effective potential, say, in lane A, 

m^Ut-Ui^ln 1 -^, (17) 

hence, the problem is essentially reduces to a one-dimensional random walk in the 
presence of the above potential field [27]. If AUi = 0, the time dependence of the 
displacement follows Sinai's law given in Eq. ( |T3l) with ip — 1/2. What makes the model 
with many particles difficult is that no effective potential exists in that case and we have 
a genuine non-equilibrium process. An interesting consequence is the following. Having 
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one particle in a finite, periodic system where the effective potential is single- valued, i.e. 
U\ = Ul + i, the expected velocity of the particle will be zero. But when putting more 
particles in the same system, the reduction to an effective potential field does not work 
anymore, and the expected velocity (or current) of particles will be non-zero (except 
when the sequence of random rates has a left-right symmetry). This interaction induced 
bias makes difficult to find a general criterion for the unbiased point J = in terms of 
the distribution of lane-change rates. Nevertheless, a sufficient condition for this is that 
the joint distribution of rates P(u, v) is symmetric under the interchange of u and v. In 
the following we will focus on this case, which has not been studied so far. 

The model has exclusively been investigated in the driven phase, where J ^ in 
finite periodic systems. It has been found that the steady state at a bottleneck with 
u > v (while in the rest of the system v > u), is much different from that of the 
one-lane partially asymmetric exclusion process PASEP. In the latter, the front of the 
high-density cluster is localized, leading to an exponentially decreasing current with the 
size I of the bottleneck, whereas in the former the front performs a symmetric random 
walk, resulting in a larger current in the order of Z _1 [27]. Based on this, the current in 
the driven phase of the disordered model was predicted to decay with the system size 
L as J(L) ~ 1/lnL, rather than the much faster algebraic decay in the corresponding 
phase of the one-lane DPASEP. Therefore, owing to the collective diffusive motion of the 
high-density cluster, particles overcome barriers more efficiently in the two-lane model 
than in the one-lane PASEP. We aim here at investigating whether this effect is present 
at the unbiased point. Since the phenomeno logical theory based on independent barriers 
breaks down here, we will resort to numerical methods. 

4-2. Steady state 

A special property of this model, which is due to a symmetry, is that if N = L, the 
stationary current is exactly zero for any set of the lane-change rates, and the local 
densities fulfill (nf) = 1 — (nf) [27]. If the global density is smaller than 1/2, the 
stationary current is non-zero (but vanishing in the limit L — > oo) and the steady state 
is segregated, consisting of a half-filled domain where approximately (nf ) + (nf ) ~ 1 
and a zero-density domain, where the density is almost zero. An appropriate measure 
of typical size of the half-filled domain is the corresponding surface width, but now the 
surface increment is related to the occupation numbers as Ahi = 2{nf + nf ) — 2N/L. As 
it is clear from the above picture, the surface width in the steady state is proportional 
to the system size, i.e. W(L,t — > oo) ~ L for iV ^ L. The steady state is, however, 
much different at half-filling (N — L). In that case, the half-filled phase covers the 
entire system, thus the mean profile (nf) + (nf ) is flat, like in a homogeneous system 
where itj = Vi — 1. As a comparison, we have investigated the homogeneous system 
by Monte Carlo simulations (data not shown) and found the surface width to obey 
the scaling relation W(L,t) = TJ-^Witj 'L 2 ) for any density, similar to the symmetric 
simple exclusion process belonging to the Edwards- Wilkinson universality class [38] . 
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Thereby one would expect the surface width in the disordered model for N = L to be 
reduced, compared to that at other densities. Surprisingly, numerical simulations still 
indicate W(L, t — >■ oo) ~ L, which refers to large fluctuations. The reason for this is the 
following. Besides the half-filled state, there are two other states, in which the current 
is trivially zero: the empty and the full lattice. In a half- filled system (N = L), let us 
imagine a state consisting of half- filled segments with {nf) = 1 — (nf ), as well as empty 
and fully occupied segments. In the homogeneous model, this state will not be stable 
(i.e. a steady state) since particles penetrate into the empty segments, just as holes into 
fully occupied segments and move there with a finite velocity for u ^ v, leading to that 
these segments ultimately dissolve and the whole system is occupied by the half-filled 
phase. In a disordered system, however, the half-filled phase is unstable. Here, fully 
occupied and empty segments are able to form spontaneously due to inhomogeneities 
and live for long times. 

To illustrate this phenomenon, let us consider a caricature of a disordered system, 
which contains just two extended regions (called reverse bias regions) of length I, where 
u > v, and which are far from each other, while in the rest of the system v > u, 
see Fig UJ An interesting property of the model under study is that a particle in an 

2 p -, 

1 1 I 1 | 

1 1 

Figure 7. Snapshot of the density profile {nf) + {nf) (green line) in a system 
consisting of two reverse bias regions (denoted by the dashed red sections). The 
arrows indicate the direction of motion of particles in the empty lattice and holes 
in the occupied lattice. 

otherwise empty region and a hole in an otherwise occupied region move in the same 
direction; from left to right (from right to left) if v > u (u > v). Let us assume, that 
a small empty segment appears around the right end point of one of the reverse bias 
regions. This segment (or condensate of holes) is stable, as particles, penetrating into 
this segment from the surrounding half-filled phase (either from left or from right) move 
preferably away from the end point, i.e. they are forced back into the half-filled phase. 
So the domain walls, separating the empty segment from the surroundings are stable. 
Of course, due to the conservation of particles, a fully occupied segment must emerge at 
the right end point of the other reverse bias region, which will be stable, as well, from 
similar reasons given above. Although the mean particle current through the half-filled 
phase in between the two reverse bias regions is zero, its fluctuations result in that the 
size of empty and fully occupied segments fluctuates (in a correlated way due to particle 
conservation), as well. The empty (or the fully occupied) region can extend over a length 
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of / on both sides of the right end point of the reverse bias region, independently. These 
correlated fluctuations of the size of particle (and hole) condensates at the two reverse 
bias regions amount to a net displacement of particles from one reverse bias region to 
the other one and lead to fluctuations of the corresponding surface. 

Concerning the inhomogeneous model, it can be divided into effective forward and 
reverse bias regions, where for the majority of links t\ > Ui and u$ < itj hold, respectively. 
Therefore, similar fluctuations occur as in the above simplified system. The largest 
effective reverse bias region is macroscopic in the unbiased case [27], hence the surface 
width must scale linearly with the system size, in accordance with the observations. 

4-3. Coarsening 

When the system is started from a homogeneous initial state, a coarsening phenomenon 
can be observed, which is analogous to that in the previous two models if the system 
is not at half-filling. In this case, the system consists of half-filled and empty (fully 
occupied) segments if N < L (N > L), the typical size of which is growing in time. 
This picture is qualitatively different for N = L, nevertheless, starting the system from 
a state with W(t — 0) = 0, the surface width must increase in time also at half-filling. 

We have performed Monte Carlo simulations for systems of size L = 2 • 10 5 , starting 
from a regular configuration where W(t = 0) = and measured the evolution of the 
surface width for not too long times, so that finite-size effects are negligible. We have 
considered different global densities N/2L = 1/2, 3/8, 1/4, 1/8, 1/16 and several samples 
for each density. The effective dynamical exponents calculated from the average surface 
width are plotted against time in Fig. |8j As can be seen in the figure, at half-filling, 




Figure 8. Effective dynamicai exponents of coarsening in the bidirectional, 
two-lane model measured in Monte Carlo simulations at different global densities 
N/2L = 1/2,3/8,1/4,1/8,1/16. Simulations have been pcrfomcd with the random 
sequential update, except one of the data sets, which was obtained by the parallel 
update for L = 48000 and 10 2 samples. 
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the effective dynamical exponent seems to tend to 0.5 but very slowly, possibly with a 
logarithmic correction. The solid curve fitting well to the data corresponds to the form 
Wit) ~ (t/lnt) 1 / 2 . Below half-filling, the effective exponent \/z e s{t) is decreasing with 
t for long times. For large enough densities, such as 3/8 and 1/4, it seems to tend to 
finite limiting values which decrease with the density. For smaller densities the effective 
exponents l/z e g(t) decrease more quickly and from the data obtained at finite time 
scales it is difficult to decide whether they tend to finite values or to zero. The latter 
case may correspond to a logarithmic coarsening law W(t) ~ [ln(i/fo)p^ with some 
finite barrier exponent ip. As a comparison, the dependence of the effective exponent 
l/z e s(t) on time in this case with ip = 0.25 and to = 1 is plotted in the figure. 

4-4- The mean- field model 

We have also studied the above model within the mean-field approximation. Here, one 
has the following evolution equations for the local densities pi(t) and vrj(t) in lane A and 
B, respectively: 



for i — 1, 2, . . . , L. The steady state of the mean- field model is, however, much different 
from that of the original one, as can be seen from the steady state of the homogeneous 
model containing a reverse bias region. The mean-field treatment cannot describe the 
diffusive motion of the half-filled particle cluster mentioned in the previous section. 
Instead, the front of this cluster is pinned in the middle of the reverse bias region. 
This leads to an exponentially decreasing current with the length of the reverse bias 
region, and to Eq. ( TT4"1) for the finite-size scaling of the current and to Eq. ( 11511 for 
the coarsening dynamics in the disordered model. We have calculated the stationary 
current at a density N/2L = 1/4 for different system sizes by solving Eqs. ( Tl8l) in 10 4 
samples for each L. The distribution of the current shown in Fig. [9] indeed follows the 
law given in Eq. (Tl~4|) . In addition to this, we have followed the evolution of the surface 
width in large systems at the same density. The time dependence of the surface width 
averaged over a few samples is again in agreement with the prediction in Eq. (I15p . as 
can be seen in Fig. [9j 

5. Discussion 

In this work, we have considered different exclusion processes, in which the disorder in 
the transition rates results in anomalous dynamics compared to the clean model. For 
the DTASEP, improving an earlier phenomenological model [12], we have obtained the 
form of logarithmic correction in the coarsening law, which is found to be compatible 
with results of Monte Carlo simulations. In the DPASEP, we have confirmed the earlier 
predicted logarithmically slow coarsening dynamics by numerical simulations. 



dpi 
dt 
diTj 

dt 



Pi-l(l - Pi) - Pi(l - p i+1 ) + ViTTi(l - Pi) - Uipiil - 71 i) 



7T i+1 (l ~ IXi) ~ 7T;(1 - TT^i) - 1^(1 - /};) + 14^(1 - 7T;) 



(18) 
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Figure 9. Left: Scaling plot of the distribution of the logarithm of the stationary 
current obtained for different sizes by the mean-field approximation. The unsealed 
data are shown in the inset. Right: Time-dependence of the surface width in the 
mean-field two-lane model. 



Then we have studied a bidirectional, two-lane exclusion process with disordered 
lane-change rates, the dynamics of which are less known than in the case of the previous 
two models. In the absence of a potential, the steady state and the non-stationary 
behavior are non-trivial and are much different compared to those of the disordered 
partially asymmetric one-lane model. In the steady state, the disorder induces formation 
of particle and hole condensates of fluctuating size, which leads to anomalously large 
fluctuations of the equivalent surface. The fact that particle clusters accumulating at 
bottlenecks are delocalized makes possible a faster coarsening than in the DPASEP. 
Numerical simulations show that for particle densities not far from 1/2, the coarsening 
is described by power-laws (apart from possible logarithmic corrections) with density- 
dependent dynamical exponents. For the particular case of the half-filled lattice, the 
dynamics are found to be characterized roughly by the exponents a = 1, (3 = 0.5 and 
z — 2, again with possible logarithmic corrections. For smaller densities, the inverse of 
the effective dynamical exponent seems to decrease constantly, but it is not possible to 
decide from the numerical data whether it tends to finite values. 

The slowing of the dynamics with decreasing density is consistent with that the 
motion of a single particle is ultra-slow, characterized formally by an infinite dynamical 
exponent. Note, however, that approaching to the zero density through positive 
densities, i.e. performing N/2L — > in an infinite system may give a limit dynamical 
law that is different from that of the single particle. 

It is worth mentioning that even for densities close to 1/2 the asymptotic region with 
a presumably finite dynamical exponent sets in only after a very long transient. In the 
transient, the coarsening, as well as the finite-size scaling of the stationary current (not 
shown) can be well described by an actived scaling with some effective barrier exponent 
xjj, that is smaller than 1/2. Probing the dynamics of the model by finite-size scaling of 
the distribution of the stationary current is hard since, due to the long transient, large 
system sizes are needed and samples in the small-current tail of the distribution have 
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Table 1. Comparison of the finite-size scaling of the stationary current in the one-lane 
DPASEP and the bidirectional, two-lane model with disordered lane-change rates. 



type of model 


one- lane 


two-lane 


single barrier of length / 
biased (F > 0) 
unbiased (F = 0) 


J ~ e~ d 
J ~ L~ a W 
— In J ~ L 1 / 2 


J~ l/l 
J~ l/ln(L) 
density dependent 



extemely long relaxation times. The mean-field version of the above model, where 
particle clusters at bottlenecks are localized, shows logarithmically slow coarsening 
formally with an infinitely large dynamical exponent. So we can conclude, that the 
slow modes related to the delocalized particle clusters near barriers are important with 
respect to the dynamics of coarsening. 

In Table HJ we have summarized the finite-size scaling behavior of the stationary 
current in disordered variants of the ASEP investigated here or earlier. Interesting 
extensions of these investigations would be the study of dynamics in multi-lane transport 
systems or in the disordered ASEP in higher dimensions [39J . 
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Appendix A. 

We have performed the Monte Carlo simulations with parallel updates according to the 
following scheme. In the case of the DTASEP, we used a parallel two-sub-lattice update 
shown in Fig. IA1I Here, the lattice is divided into even and odd sub-lattices, which are 
updated alternately. In the case of the DPASEP, a similar two-sub-lattice update has 
been applied as given above. 



Parallel threads: 1,2,.. 
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Figure Al. Illustration of the parallel, two-sub-lattice update applied for the 
DTASEP 



Anomalous coarsening in disordered exclusion processes 



16 



In the case of the bidirectional, two-lane model with parallel update each Monte 
Carlo step consists of the following operations: (i) even sub-lattice updates of the lanes 
A and B, with probabilities 1/2, (ii) lane changes with probabilities Ui and Vi, (iii) odd 
sub- lattice updates of the lanes A and B with probabilities 1/2, (iv) lane changes with 
probabilities Ui and Vi. Further details can be found in Ref. [23] . 
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